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ABSTRACT 


The  objective  of  this  work  is  to  study  the  transfer  of  energy 
between  an  elastic  solid  coating,  a  viscoelastic  non-homogeneous 
layer  and  the  surrounding  viscous  fluid.  The  particle  displacements 
and  stress  distributions  in  the  solid  and  the  fluid  are  analyzed.  The 
motion  at  the  lower  boundary  of  the  elastic  solid  is  assumed  given. 
The  effect  of  coating  properties,  transition  layer  thickness  and 
frequency  are  obtained. 


I  -  INTRODUCTION 


The  drag  of  a  flow  field  over  a  surface  will  be  reduced  if 
the  energy  content  of, the  fluid  can  be  decreased.  The  design  of 
a  drag  reducing  compliant  coating  thus  requires  a  mechanism  where¬ 
by  the  coating  interacts  with  the  fluid  in  such  manner  as  to  either 
absorb  or  cancel  the  fluid  disturbances.  The  accomplishment  of 
this  goal  requires  an  understanding  of  the  transfer  mechanism  be¬ 
tween  the  solid  and  the  fluid. 

The  transfer  of  energy  between  a  solid  and  a  fluid  is  accom¬ 
plished  by  the  continuity  of  particle  displacement  and  stresses 
across  the  interface-*-.  The  presence  of  a  viscoelastic  non-homo- 
geneous  layer  interposed  between  a  solid  and  a  fluid  (Fig.  1)  with 
properties  varying  from  that  of  the  solid  to  that  of  the  fluid  at 
the  two  extremes  obscures  the  actual  identity  of  the  interface  and 
of  the  boundary  conditions. 

Many  different  reasons  exist  for  justifying  the  interest  in 
the  analysis  of  a  thin  viscoelastic  non-homogeneous  transition 
layer  between  a  coating  and  the  surrounding  fluid.  Kramer2  in 
some  of  the  early  compliant  coating  work  applied  a  thin  layer  of 
paint  on  the  rubber  surface  of  the  coating.  From  a  physical  point 
of  view,  thi«  transition  layer  can  "bridge  the  gap"  between  the 
solid  and  the  fluid  in  a  more  gradual  manner  and  this  may  alter 
the  energy  transfer  between  the  two  media.  When  a  compliant  sur¬ 
face  is  exposed  to  a  turbulent  flow  field,  the  surface  will  experi¬ 
ence  forces  which  will  vary  randomly  in  both  the  space  and  the  time 
domain.  The  resulting  surf  ace  motion  will  be  composed  of  random  modes 
superimposed  upon  some  form  of  discrete  motion.  This  complex 
motion  will  cause  the  classical  wall  boundary  condition  to  be 
applied  at  a  mean  surface  location  which  may  vary  as  a  function 
of  both  time  and  space.  The  set  of  interface  boundary  conditions 
may  then  be  modelled  by  a  system  with  the  classical  wall  boundary 
conditions  being  satisfied  along  a  layer  whose  thickness  is  equiva¬ 
lent  to  the  crest  to  trough  distance  of  the  mean  motion  (Fig.  1) . 
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Two  different  modes  of  energy  transfer  between  a  compliant 
coating  and  the  surrounding  fluid  exist:  shear  and  longitudinal 
waves.  This  report  restricts  its  attention  to  shear  waves  only. 
For  the  sake  of  simplicity,  the  motion  will  be  assumed  to  be  one¬ 
dimensional  (Fig.  2)  and  wall  conditions  prescribed  at  the  lower 
boundary  of  the  solid. 


II.  ANALYSIS 


II-l  Derivation  of  Governing  Equations 

Consider  a  plane  strain  problem  with  no  variations  in  the  2 
direction  (Fig.  2) .  The  general  equations  of  motion  valid  for  a 
general  medium  are3: 
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where  v.  and  v,  are  the  velocities  in  the  x  and  the  y  direction 
^  ^  D 

respectively  and  ^  is  the  material  derivative  given  by 
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For  a  one-dimensional  (no  variations  in  the  x  direction)  steady 
state  (exp(-icut))  system,  the  governing  equations  simplify  to: 
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If  only  shear  modes  are  considered,  then  v^  -  0 ;  the  system  of 
differential  equations  is  reduced  to  one. 
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+  ip  a:  v^  =0 


(18) 


where  x-j^  is  9iv©n  by  either  Eq.  (14)  or  Eq.  (15).  The  stress-strain 
relationships  simplify  for  a  fluid  with  zero  mean  flow  velocity  and 
with  steady-state  disturbances  given  by 


X1  =  U1  exP(~ia)t) 


(19) 


The  stress-strain  relationships  for  the  solid  and  the  fluid  become 


Fig.  3  shows  the  magnitude  of  the  stress-strain  coefficient 

as  a  function  of  axial  distance.  In  a  typical  solid-fluid  system, 

the  magnitude  of  the  coefficient  jumps  from  a  value  of  G  at  the 

solid  to  a  magnitude  of  uw  in  the  fluid.  The  lowest  values  of  the 

4  2 

shear  modulus  for  a  very  soft  plastisol  is  about  10  N/m  .  For 

4  2 

water  at  10  Hz  the  value  of  uoj  is  about  60  N/m  .  For  more  rigid 

surfaces  and  for  lower  frequencies,  the  difference  between  the 

value  of  G  and  uw  increases. 

In  the  presence  of  a  transition  layer  the  change  can  occur  over 
a  finite  thickness.  Introducing  a  generalized  stress-strain  relation¬ 
ship  will  allow  one  to  combine  ”qs .  (20)  and  (21)  into  a  single  equa¬ 
tion.  This  new  unified  approach  leads  to 

dul 
dy 

where  K  is  equal  to  G  within  the  coating  material,  and  is  equal  to 
-iuu)  in  the  fluid;  within  the  transition  layer  K  takes  on  intermediate 
complex  values.  The  governing  differential  equation  then  reduces  to: 


=  K 


(22) 


h  [K  ly>  sr] '+  °“2ui  * 0  •  1231 

For  the  particular  case  of  interest  in  this  study,  where  the 
interface  region  between  the  elastic  solid  and  the  fluid  is  assumed 
to  be  a  non-homogeneous  viscoelastic  layer,  the  value  of  K (y)  is 
assumed  to  be  given  by 

K_(y)  =  — i +  £g  +  iyuj  exp(-y/Z)n  ,  (24) 

where  l  is  the  thickness  of  the  compliant  coating  and  n  is  a  parameter 
inversely  proportional  to  the  thickness  of  the  transition  layer.  Note 
that  for  y <<Z  one  recovers  the  shear  modulus  of  the  solid  while  for  y»Z 
one  obtains  the  fluid  viscosity  relationship.  For  simplicity  it  is 
assumed  that  the  solid  and  the  fluid  have  equal  densities .  The  value 
of  K (y )  is  non-dimens ionali zed  with  respect  to  G,  thus 

3y  (H  3y^)+  (r2p)  U1  "  0  '  (25) 

where 

H(y)  -  -iR  +  (1  +  iR)  exp(-y/Z)n  .  (26) 

The  value  of  R  is  defined  to  be  the  ratio  of  uus  to  G.  The  shear 
wave  speed  in  a  solid  is  given  by 
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c2s  =  G/p  ,  (27) 

while  in  a  fluid  the  wave  speed  of  the  shear  waves  is  given  by 
2 

cf  =  u<i)/p  (28) 

The  value  of  P  can  then  be  written  as  the  ratio  of  the  two 
speeds 

R  =(cf/cs)2  =  Uf/Xg)2  -  (29) 

Fig.  4  shows  the  variation  of  the  magnitude  of  H  as  a  function 
of  y/i  for  different  values  of  the  layer  thickness  parameter  n.  Note 
that,  in  the  limit  as  n-1-®  ,  the  classical  no-transition  layer  config¬ 
uration  is  recovered. 

The  values  of  distance  y  and  displacement  u^  have  been  non- 
dimensionalized  with  respect  to  the  coating  thickness  l  and  the 
magnitude  of  wall  displacement  respectively.  Eq.  (25)  together  with 
the  proper  boundary  conditions  represent  the  equation  to  be  studied 
in  this  report.  The  boundary  conditions  of  the  problem  are  split; 
one  condition  is  at  the  lower  boundary  of  the  coating  while  the  t 

other  is  the  radiation  condition  at  large  distances  from  the  sur¬ 
faces.  The  inner  boundary  of  the  coating  can  move  by  imposing  either 
a  displacement  or  a  force.  The  split  boundary  conditions  then  become 

u1(0)  *  1  , 

lim  u^y)  =  C  expf-d+i)  ry/(/?i)l  ,  (30) 

y-4>oo  ^  ^ 


-  1  / 
y~o 

u1(y)  =  C  exp  [-(1+i)  ry/(/2£)]  .  (31) 

Boundary  conditions  as  represented  by  Eq.  (30)  will  give  dis¬ 
placement-driven  solutions  while  conditions  of  Eq.  (31)  will  give 
force-driven  solutions.  Note  that  since  the  stress-strain  coeffic¬ 
ients  have  been  non-dimensionalized  with  respect  to  G,  the  stress 
at  the  lower  boundary  is  equal  to  the  derivative  of  the  particle 
displacement. 


UU^ 

W 

lim 

y-M» 


7 


II-2  Numerical  Solution. 

The  solution  to  the  problem  considered  in  this  report  requires 

the  numerical  integration  of  Eq.  (25) 
o  du,  - 

^  (H  +  (I^R)  u1  =  0  .  (25-a) 

The  boundary  conditions  for  the  displacement-driven  case  are 
Eq.  (30) 

u, (0)  =  1 

lim  ux(y)  =  C  exp  [-(1+i)  Fy/ (/?£)]  .  (30-a) 


Eq.  (25-a)  was  solved  numerically  by  the  Runge-Kutta  technique 
Values  of  slopes  of  displacements  at  y=0  were  iterated  until  the 
far-field  boundary  condition,  Eq.  (30-a) ,  was  satisfied.  For  most  . 
cases  of  interest  this  procedure  required  5  t:o  7  iterations.  The 
iterations  were  performed  by  studying  the  sensitivity  of  the 
assumed  slope  at  the  wall  (note  that  the  values  of  the  displace¬ 
ments  and  slopes  are  complex)  upon  the  far-field  radiation  condition 
The  -values  of  the  parameters  rand  R  can  be  re-written  in  terms 
of  length  ratios  as 


r  =  2  it  u/\f) 

r  =  <vv2 

The  range  of  values  of  T  and  R 
immersed  in  water  are  shown  in 


(33) 

for  different  coating  materials 
Table  1. 


Cg (m/sec) 

R 

2500 

IQ'12 

250 

lo-io 

25 

10~8f 

2.5 

10  ~6f 

I  (mm)  T 

100  250  /F 

10  25  /f 

1  2.5/F 

0.1  0.25/f 


Table  1-  Values  T  and  R  for  different  coatings  • 


s 


II-3  Modified  Stokes  Second  Problem  (MSS?) 

In  the  limit  as  Hoth  the  transition  layer  and  the 
elastic  solid  layer  thicknesses  go  to  zero  (i.e.  n-*"»  and  2.-0 )  ,  one 
will  obtain  the  classical  Stokes  second  problem**.  Stokes  second 
problem  corresponds  to  the  propagation  of  shear  disturbances  in  a 
viscous  fluid  produced  by  an  oscillating  plane  wall.  For  a  zero¬ 
thickness  transition  (n-*-°°)  but  a  finite-thickness  elastic  solid  one 
obtains  the  modified  Stokes  second  problem  (MSSP) .  In  the  MSSP, 
instead  of  driving  the  viscous  fluid  with  a  given  displacement 
oscillation,  the  lower  boundary  of  the  solid  is  forced  by  the  same 
function.  The  displacement  on  the  lower  boundary  causes  shear  waves 
to  be  transmitted  first  through  the  elastic  medium,  then  into  the 
fluid.  The  fluid-structure  coupling  occurs  at  the  solid-fluid  inter¬ 
face  where  continuity  of  displacement  and  stresses  are  required. 

The  equations  of  motion  representing  the  MSSP  are 

32u.  3t. 

P  =—  =  y-—  0  <v  <  Z  ,  (34) 


where  the  solid  and  the  fluid  stresses  are  related  to  the  particle 
displacement  u^  by 


T1  =  G 


°iyi£ 


T,  -  U  ItT  y-£  1  (37) 

1  oy  3 1 

The  boundary  and  the  interface  boundary  conditions  are  given 


ul<0> 

*  1  exp(-iu)t) 

U.  U) 

j  =  u,  (£)  1 

X 

1  solid  A  1 

fluid 

T,  (£) 

1  =  T. (t)  j 

X 

'solid  1  1 

fluid 
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The  solution  to  Eqs .  34-35  subject  to  boundary  conditions 
38-40  can  be  obtained  in  closed  form.  It  can  be  shown  that  the 
values  of  the  particle  displacement  un  in  the  two  regions  is  given 


ux(y) 
ui  (y) 


it  j, 

[cos  S  -  (R/2)  ~  sin  3  ]-  i  (R/2)  ~  sin  ? 

D 

exp[  (i-1)  r(y-n /(£/?)  1/D  y>* 


0  <y<Z  (41) 


(42) 


where 

D  =  cos  ft  -  (R/2)  55  sin  9  -  i  (R/2)*  sin  ft,  (43) 

3  =  ©  (1-y /l)  ,  (44) 

9  =  T/R  =  2tt£  /  A  .  (45) 

s 

Note  that  in  the  limit  as  ©-*-0  one  recovers  the  solution  to 
the  classical  Stokes  second  problem  given  by 

uL(y)  *  exp[  (i-1)  T  (y-2.)  /  (1/2  )]  .  (46) 

Two  other  classical  solutions  are  recoverable  from  Eqs.  (41) 
and  (42).  In  the  case  of  an  inviscid  fluid  (i.e.  R-*-0)  the  solution 
reduces  to 

uL(y)  =  cos 8/cos  0  .  (47) 


In  the  case  of  a  very  viscous  fluid  (i.e.  R-*°°)  the  solution 
reduces  to 

u^  =  sin  0/sin  9  .  (48) 

Eqs.  (47)  and  (48)  correspond  to  the  solutions  for  shear  waves  in  a 
finite- thickness  elastic  medium  with  a  free  and  a  fixed  boundary. 
Conditions  for  9  =  (l+2n)tr/2  and  9  =  2n~  for  the  former  and  the 
latter  case  correspond  to  the  resonance  conditions  associated  with 
a  closed-open  and  a  closed-closed  pipe  respectively  ^ . 

The  solution  to  the  MSSP  as  represented  by  Eqs.  (41)  and  (42) 
together  with  the  limiting  cases  (4,7)  and  (4  8)  will  become  very 
useful  in  understanding  the  results  obtained  in  the  presence  of  a 
transition  layer. 
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III  -  DISCUSSION  OF  RESULTS 


III-l  Displacement -Driven  Boundary 

In  the  discussion  that  follows,  the  results  are  presented  in 
the  form  of  particle  displacement  and  shear  stress  distribution  as 
function  of  distance  from  the  lower  boundary;  a  unit  wall  distur¬ 
bance  is  applied  at  the  lower  boundary  of  the  coating.  The  numer¬ 
ical  integration  of  the  equation  of  motion  for  the  conditions  of 
very  thin  transition  layer  (i.e.  n-*-«°)  should  give  results  which 
agree  with  the  solution  to  the  MSSP  discussed  in  the  previous 
section.  Figures  6  and  7  show  the  particle  and  the  shear  stress 
distribution  for  typical  conditions  of  four  increasingly  thinner 
transition  layers.  The  numerical  solutions  are  seen  to  asymptote 
to  the  proper  limiting  solution.  The  transition  layer  for  the  four 
values  of  n  considered  (n  =  8,  16,  32,  64)  correspond  to  thicknesses 
of  about  0.54  ,  0.254  ,  0.124,  and  0.064  respectively.  It  is  inter¬ 
esting  to  note  that  as  the  transition  layer  thickness  increases,  a 
particle  displacement  overshoot  (i.e.  local  displacement  larger 
than  unity)  begins  to  appear.  The  overshoot  becomes  more  evident 
as  either  T  or  R  is  increased.  Physically  these  conditions  corres¬ 
pond  to  higher  frequency ,  thicker  elastic  layer  and  softer  coating 
material  (lower  Cg) .  Figs.  8  and  9  show  the  effect  of  thicker 
elastic  layers  (bigger  r  )  while  Figs.  10  and  11  show  the  effect 
of  softer  coating  materials  (higher  R) .  In  both  instances  the 
transition  layer  thickness  has  been  fixed  to  correspond  to  n  =  8. 
Reducing  the  transition  layer  thickness  will  cut  down  on  the  magni¬ 
tude  of  the  overshoot  but  will  be  eliminated  only  in  the  limit  as 
n-*-°°.  A  transition  layer  equal  to  0.54  (n=8)  will  produce  a  48% 
overshoot  while  a  0.254  thick  layer  will  produce  a  20%  overshoot. 
These  values  are  for  R*  10~5  and  r  =*  50  (Fig.  12)  .  The  no -transit  ion  layer 
results  agree  exactly  with  the  present  analysis  at  all  locations 
except  for  values  of  y  in  the  neighborhood  of  the  overshoot.  The 
overshoot  produced  by  the  classical  approach  (i.e.  no  transition 
layer)  is  minimal  and  it  cannot  be  seen  on  the  scales  used  in  the 
figures.  The  relative  magnitude  of  the  overshoot  is  a  function  of 
the  transition  layer  thickness  and  elastic  layer  properties. 


The  physical  reasoning  for  the  overshoot  can  be  best  understood 
by  referring  to  the  analytical  solution  of  the  MSSP  (Zqs.  41  and  42). 
Evaluating  the  particle  displacement  at  the  solid-fluid  interface  gives 


U  x(£)  =  D_1  .  (49) 

Values  of  T  and  R  which  will  minimize  D  will  maximize  the 
particle  displacement  at  the  interface  and  in  its  neighborhood. 

It  is  instructive  to  note  that  in  the  limiting  condition  of  an 
inviscid  fluid  (i.e.  R-*-0)  and  with  9  =  it/2  (i.e.  coating  thickness 
1  equal  to  one  quarter  of  shear  wavelength  of  solid)  one  will 
obtain  very  large  particle  displacements  at  the  interface.  This 
result  is  consistent  with  the  concept  of  a  velocity  (or  displacement) 
driven  open  pipe  having  waves  reflected  at  the  open  end  of  the  pipe 
and  creating  resonance  conditions.  The  displacement  overshoots  are 
thus  seen  to  be  caused  by  the  reflection  of  the  shear  waves  at  the 
solid-fluid  interface.  The  presence  of  a  viscoelastic  non-homo- 
geneous  transition  layer  at  the  interface  is  seen  to  magnify  the 
internal  reflection  of  the  waves.  The  actual  magnitude  of  the 
magnification  is  a  function  of  the  three  parameters  T,  R  and  n. 

III-2  Force-Driven 

For  conditions  where  the  lower  boundary  of  the  elastic  layer 
is  force-driven,  the  wall  boundary  conditions  applicable  to  the 
numerical  integration  of  Eq.  25  is 

-  1  .  ( 31-a) 

y=0 

The  value  of  the  wall  displacement  u^O)  is  iterated  until  the 
proper  far-field  radiation  condition  is  obtained. 


\ 


Two  typical  cases  of  force -driven  conditions  are  shown  in  Figs. 
13-15.  When  the  displacement  profiles  are  non-dimensionalized  with 
respect  to  the  wall  displacements,  the  results  collapse  on  the  same 
curve  obtained  for  the  displacement-driven  condition  (see  Fig.  16)  . 
The  same  conclusion  applies  to  the  shear  stress  distribution.  It 
should  also  be  noted  that,  as  expected,  the  Rayleigh  reciprocal 
theorem  =  T2U]_)  holds.  This  result  can  be  verified  by  noting 

that  at  any  y  location  the  product  of  the  stress  from  the  dislace- 
ment-driven  case  multiplied  by  the  displacement  from  the  force  - 
driven  condition  equals  the  product  of  the  stress  from  the  force- 
driven  case  with  the  displacement  from  the  displacement-driven 
condition. 


IV  -  CONCLUSIONS 


The  effect  of  a  viscoelastic  non-homogeneous  transition  layer 
located  at  the  interface  between  an  elastic  solid  and  a  viscous 
fluid  has  been  analyzed.  The  displacement  and  the  stress  distribu¬ 
tion  due  to  either  a  displacement  or  a  stress-driven  lower  boundary 
have  been  analyzed.  It  has  been  shown  that  a  thin  layer  can  radically 
change  the  reflection  characteristics  of  the  interface  in  such  manner 
as  to  magnify  the  displacement  overshoot  that  is  normally  present. 

Results  have  been  presented  and  discussed  for  displacement  and 
for  force-driven  coating  disturbances.  No  basic  differences  exist 
between  the  two  solutions.  The  present  results  have  been  analyzed 
with  and  compared  to  the  analytic  solution  of  the  modified  stokes 
second  problem. 
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FIG.  i  -  TRANSITION  LAYER  MODELS  -  (a)  viscoelastic 

NON* HOMOGENEOUS  LAYER  -  (b)  EQUIVALENT  TRANSITION 
LAYER  DUE  TO  SURFACE  MOTION. 
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FIG.  E  -  ONE-DIMENSIONAL -TRANSVERSE  WAVES  IN 
SOLID;  TRANSITION  LAYER,  AND  FLUID. 
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FIG.  3-  VARIATION  OF  THE  STRESS- STRAIN  COEFFICIENT 
AS  A  FUNCTION  OF  WALL  DISTANCE  FOR  (a)  NO 
TRANSITION  ,  (b)  WITH  TRANSITION  LAYER. 
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4  -  TRANSITION  LAYER  CHARACTERISTICS.  R  =  IQ3 


PARTICLE  DISPLACEMENT  DISTRIBUTION  -  EFFEC 
OF  TRANSITION  LAYER  THICKNESS.  R=IO'5,r  =  . 

- - OUTER  EDGE  OF  TRANSITION  LAYER. 


SHEAR  STRESS  DISTRIBUTION  AS  FUNCTION  OF  TRANSITION 
LAYER  THICKNESS.  R=IO'sr  =  3. 


thicker 

COATING 


FIG.8-PARTICLE  DISPLACEMENT  DISTRIBUTION 

FOR  DIFFERENT  COATING  THICKNESSES 
R=IO"5,  r\«8. 


STRESS  DISTRIBUTION  AS  FUNCTION  OF  COATING  THICKNESS. 
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FIG.  10-PARTICLE  DISPLACEMENT  DISTRIBUTION  -  EFFECT  OF  COATING 
PROPERTIES.  n=  8,  T  =  3  . 
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FIG.  13”  PARTICLE  DISPLACEMENT  DISTRIBUTION  FOR 
FORCE-  DRIVEN  BOUNDARY.  r\=8)R=!0"5,  T=I0. 
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FIG.  14-  PARTICLE  DISPLACEMENT  DISTRIBUTION  FOR 
FORCE-DRIVEN  BOUNDARY,  a* 8  ,  R=  I0"5,  P=50. 
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FIG.  16-  NORMALIZED  PARTICLE  DISPLACEMENT 

DISTRIBUTION  FOR  FORCE- DRIVEN  BOUNDARY. 
n  =  8,  R-IO'5. 
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APPENDIX 


Enclosed  are  additional  figures  of  particle  displacement  and 
shear  stress  distribution. 
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F1G.A-I  -  PARTICLE  DISPLACEMENT  DISTRIBUTION 

FOR  DIFFERENT  COATING  THICKNESSES. 
n-8,  R  =  IO"7 . 


,.A-a- SHEAR  STRESS  DISTRIBUTION  AS  FUNCTION  OF  COATING  THICKNESS. 
n=8,  R=  10”7. 


FIG. A-3  -  PARTICLE  DISPLACEMENT  DISTRIBUTION  AS 
FUNCTION  OF  FREQUENCY  FOR  TYPICAL 
COMPLIANT  COATING.  n=8J  1.25  mm ,  Cs=0.87sec 


FIG.A-4  -  PARTICLE  DISPLACEMENT  DISTRIBUTION  AS 
FUNCTION  OF  FREQUENCY  FOR  TYPICAL 
COMPLIANT  COATING.  n«  8,  i  =  2.5mm, c,=  0.1  m/se 


SHEfcR  STRESS  DISTRIBUTION  AS  FUNCTION  OF  FREQUENCY  FOR 

TYPICAL  COMPLIANT  COATING.  n=8, 1=  2.5  mm ,  Cs«  0.1  m/sec  . 
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